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We study the kinetics of spin coherence of optically excited electrons in an undoped insulating 
ZnSe/Zni-zCd^Se quantum well under moderate magnetic fields in the Voigt configuration. After 
clarifying the optical coherence and the spin coherence, we build the kinetic Bloch equations and 
calculate dephasing and relaxation kinetics of laser pulse excited plasma due to statically screened 
Coulomb scattering and electron hole spin exchange. We find that the Coulomb scattering can 
not cause the spin dephasing, and that the electron-hole spin exchange is the main mechanism of 
the spin decoherence. Moreover the beat frequency in the Faraday rotation angle is determined 
mainly by the Zeeman splitting, red shifted by the Coulomb scattering and the electron hole spin 
exchange. Our numerical results are in agreement with experiment findings. A possible scenario for 
the contribution of electron-hole spin exchange to the spin dephasing of the n-doped material is also 
proposed. 
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I. INTRODUCTION 

Studies of ultrafast nonlinear optical spectroscopy in 
semiconductors have attracted numerous interest both 
experimentally and theoretically during the past 20 
yearsBH Most of these studies are focused on the 
optical coherence and the studies of spin coherence 
are relatively rare. _ . Recently, ultrafast nonlinear op- 
tical experimentsatLia have shown that the spin co- 
herence which is optically excited by laser pulse can 
last much longer than optical coherence. For undoped 
ZnSe/ZnCdSe quantum wells, it is found in the experk 
ment that the spin coherence can last up to 15-20 psci 
where as for undoped bulk GaAs, it lasts about 600 ps.t2 
For n-doped material, the spin coherence can last up to 
three orders of magnitude longer than the undoped sam=. 
pie which makes 8 ns for .ZjiSe/ZnCdSe quantum welltll 
and 100 ns for bulk GaAs.liia These discoveries have given 
rise to an emerging interest within the physics and elec- 
tronic engineering communities in using electronic spins 
for the storage of coherence and also have stimulated the 
optimism that the coherent electrons will finally be real- 
ized as a basis for quantum computation. 

The electron spin coherence can be directly observed 
by femtosecond time-resolved Faraday rotation (FR) in 
the Voigt configuration. In that configuration a moderate 
magnetic field is applied normal to the growth axis of the 
sample. The coherence is introduced by a circularly po- 
larized pump pulse that creates electrons and holes with 
an initial spin orientation normal to the magnetic field. 
Then the magnetic field causes the electron spin to flip 
back and forth along the growth axis which makes the 
net spin precess about the magnetic field. The hole spin 
is kept along the growth axis direction of the quantum 
well as will be discussed below. After a certain delay 
time t, a linearly polarized probe pulse is sent into the 
sample along a slightly different direction from the pump 
pulse and by measuring the FR angle, one can sensitively 



detect the net spin of electrons associated with the delay 
time t. This method has proven to be extremely suc- 
cessful in measuring the coherent spin evolution and spin 
dephasing»Ery 

In order to further extend the spin coherence time, it 
is important to understand the physics of spin dephas- 
ing. While there is extensive theoretical study and un- 
derstanding of the optical dephasingjiil the theoretical 
investigation on the spin dephasing-is relatively limited, 
nevertheless of much longer history.L3 The early work in- 
cludes that, of Elliott in 1954, who discussed the spin 
relaxation!!^ induced by lattice and impurity-scatterings 
by taking into account the spin-obit effectsHI Later, in 
1975 Bir et al. calculated the spin dephasing using a 
model Hamiltonian describing the electron-hole spin ex- 
change (EHSE) by considering Coulomb scattering be- 
tween electron and hole, combined with the spin-obit- 
coupling induced band mixingJiB In that paper by using 
Fermi Golden rule, Bir et al. proposed a spin relaxation 
rate which is proportional to the hole density. It was 
not until 90 's that experimentalists found such an effect 
and claim the EHSE is.important in the intrinsic and p- 
doped semiconductors For n-doped samples, however, 
as the density of the electrons is much higher than that 
of the holes, the holes recombinc with electrons in a time 
much shorter than the measured spin dephasing time. 
As the predicted spin depliasipg rate induced by EHSE 
is proportional to hole density,ta it is therefore suggested 
that foij-^n-doped sample, the dephasing mechanism is 
unclearO Recently Linder and SharrO presented a the- 
ory of the spin coherence of excitons by studying Bloch 
equations. However, they did not discuss the dephas- 
ing mechanisms explicitly, and instead described all the 
dephasing by using phenomenological relaxation times. 

In this paper, we present a model to study the kinet- 
ics of spin precession of a femtosecond laser-pulse ex- 
cited dense plasma in a quantum well in the framework 
of the semiconductor Bloch equations combined with 
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carrier-carrier scattering in the Markovian limit. Non- 
Markovian effects are not important in our case, as the 
width of the laser pulse is large (100 fs) and the time scale 
of spin dephasing is very long. The main purpose of this 
paper is to understand the spin decoherence. It has been 
well known that both carrier-carrier Coulomb scattering 
and carrier-phonon scattering play significant roles in the 
optical dephasing. For the|2D carrier density around 10 11 
cm~ 2 in the experiment, til the main carriers are electron- 
hole plasma and the Coulomb scattering gives a fast de- 
phasing of the optical coherence with dephasing times, 
ranging from tens of femtoseconds to subpicoseconds,cl 
depending on the strength and width of the pump pulse 
and/or the density poLdoping, which affects the building 
up of the screening .oEJ Besides the fast dephasing due to 
Coulomb scattering, carrier-phonon scattering also con- 
tributes to the optical dephasing with the dephasing time 
being around ten picoseconds and carrier-density inde- 
pendent. For spin coherence, the experimental evidence 
that the snia_dephasing depends strongly on the car- 
rier densitjO'Ej clearly rules out the possibility that the 
main mechanism of the spin dephasing is due to carrier- 
phonon scattering. Therefore we focus on the effect of 
carrier-carrier scattering. We distinguish the spin co- 
herence from the optical coherence and study the roles 
of Coulomb scattering and EHSE to the spin dephasing 
separately. We find that the pure Coulomb scattering — 
although it destroys the optical coherence strongly — does 
not contribute to the spin dephasing at all, and that 
EHSE is the main mechanism of the spin decoherence. 

Our paper is organized as follows: we present our 
model and kinetic equations in Sec. II. Then in Sec. 
Ill we present the numerical results for an undoped 
ZnSe/Zn 1 _ :r Cd ;c Se quantum well. A conclusion of our 
main results and a discussion of the theory for the n- 
doped sample are given in Sec. IV. 



II. MODEL AND KINETIC EQUATIONS 



A. Model and Hamiltonian 



We start our investigation of a quantum well with its 
growth axis in the z-direction. A moderate magnetic field 
B is applied in the rr-direction. Landau quantization is 
unimportant for the magnetic field B in our investiga- 
tion. We consider the conduction band (CB) and the 
heavy hole (hh) valence band (VB). Due to the presence 
of the magnetic field, the spins of electrons and holes are 
no longer degenerate and therefore each band is further 
separated into two spin bands with spin ±1/2 for elec- 
trons in the CB and ±3/2 for those in the hh VB. It is 
noted that these spin eigenstates are defined with respect 
to the ^-direction. 

In the presence of the moderate magnetic field and 
with the interactions with a coherent classical light field, 



the Hamiltonian for electrons in the CB and VB is given 
by 



H = Y1 £ ^ c \ka C ^ka + • ^ S ^a'cl kr7 C^ lko 

iika 



aa' 



He + H i , 



(i) 



with fi = c and v standing for the CB and the VB 
respectively. e vk is the energy spectrum of an elec- 
tron in the VB (CB) with momentum k. It is noted 
that k here stands for a momentum vector in the x-y 
plane. e vk = ~E g /2 - k 2 /2m h = -E g /2 - e hk and 
e c k = Eg/2 + k 2 /2m e = E g /2 + e e k with rrih and m e 
denoting effective masses of hh and electron separately. 
Eg is unrenormalized band gap and a is the spin index. 
For electron in the CB, a — ±1/2 and for electron in 
the hh VB, a — ±3/2. [Ib is Bohr magneton. S M are 
the spin matrices with S c being spin 1/2 matrices for 
electrons and S v being spin 3/2 matrices for holes. 

H E in Eq. (|l]) denotes the dipole coupling with the 
light field E a (t) with a = ± representing the circular 
polarized light. Due to the selection rule the electrons 
in the spin 3/2 (—3/2) hh band can only absorb a left 
(right) circular polarized photon and go to the spin 1/2 
(-1/2) CB. Therefore 

H E = -dJ2[E-(t)cl ki c vki + H.c] 

k 2 

-d^E+ityl^^.+H.c] . (2) 

k 2 

In this equation d denotes the optical-dipole matrix el- 
ement. The light field is further split into E a (t) = 
E°(t) cos(wt) with u> being the central frequency of the 
coherent light pulse. E®(t) describes a Gaussian pulse 
E®e~ l / St with St denoting the pulse width. 

Hi is the interaction Hamiltonian. As said before, we 
focus on the carrier-carrier scattering. Therefore Hj is 
composed of Coulomb scattering and EHSIiO with thn 
later being much weaker than the Coulomb scattering.tJ 



Hi o ^g C lfc+oer C lfc' -ga' c vk'a' c ^ka 



aa' 



n UqCTa C^ k+qa cl k , _ qa ,C„k'(j' 



aa' 



lY. U 'l <7a ' C U+ q * C lk>- q a> C Hk>aC„ka> ■ (3) 



aa' 



The first term of Eq. (J3|) is the ordinary Coulomb inter- 
action. The second term describes the "direct" EHSE 
which scatters an electron in the band \x with spin er and 
momentum k to the same band with spin a and momen- 
tum k+q and in the moan time which scatters an electron 
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from the different band v /i) with spin a' and momen- 
tum k! back to that band with spin a' and momentum 
k' — q. The last term is "exchange" EHSE which scatters 
an electron in the band p with spin a and momentum 
k' to the different band v (=/= p) with spin a' and mo- 
mentum k' — q and in the mean time which scatters an 
electron from the band v with spin a' and momentum k 
back to band p with spin a and momentum k + q. It is 
noted here that a and a' stand for ±1/2 when they are 
at the CB and ±3/2 when they are at the VB. As the 
exchange effect involves a form factor which consists the 
overlap of the wavefunctions of the CB and VB, U' q is 
much smaller than U q . Moreover, it will be shown later 
that the exchange effect is also energetically unfavorable 
as it scatters carriers across the band. Therefore the di- 
rect EHSE is the dominant effect. 

For Voigt configuration, as B is along the x direction 
and (S^)±|.±| = 0, one can see directly from Eq. (|l|) 
that the spin of hh cannot be flipped by the magnetic 
field. We point here that this is only true when the mag- 
netic field is small compared to the hh-light hole split- 
ting and hence the flip between the hh and the light hole 
can be neglected in the timescale of the discussion. As 
ga; _ a x J 2 w jth a x standing for Pauli matrix, it is there- 
fore straightforward to see from the Hamiltonian that the 
magnetic field causes the CB electron spin to flip and flop. 

It can be seen directly from the Hamiltonian He that 
the laser pulse introduces the optical coherences into the 
system which are built between the CB and the VB with 
same spin direction: Vk\\ = ( c j^3 c cfci) an d Vk-\-\ = 
(c^ fc _ 3 c ck _ i ). In the mean time, due to the presence 
of the magnetic field in Eq. (|l|) , these optical coherences 
may further transfer coherence to p k _xi = ( c ^ fc 3 c c fc-i)j 
p fe i_| = {c[ k _^c ck i) and p ccfc i_i = (c| t _ic cl i>, with 

the first two being the coherence between the CB and 
VB with opposite spin directions and the last one being 
the coherence between two CB's with opposite spin di- 
rections. While it is well known that optical coherence is 
represented by p k i 3 and p k _ i _ 3 , we will show later that 
the spin coherence of electrons is represented by p cck i _ i . 
When there is no dephasing effect added on this term, 
the electron spin precession will last forever. Moreover, 
when p cck i_i decays to zero, there is no electron spin 
precession. Therefore we refer it as spin coherence in 
the following. Pu 1 3 and P k i 3 , which describe the 
coherence between the states with the optical transition 
between them being forbidden by the selection rule, may 
play certain role in the optical dephasing and we will refer 
them hereafter as forbidden optical coherences. Finally 
there is also a coherence between two VB's with opposite 
spin directions: p vvk s_^ = (c ,_3C t u}. In the absence 

2 2 VK 2 2 

of spin flip of the hh, this coherence is much weaker than 
the other coherences as it then can only be excited by 
the laser pulse through the coupling to the forbidden co- 
herence. Therefore in the present discussion of undoped 



material, we do not include it in our model. 



B. Kinetic equations 

We build the semiconductor Bloch equations for 
the quantum well by nonequilibrium Green function 
methodE-3 as follows: 



PfJ,v,k,acr' — P[iv,k,aa' Icoh P^v.k.aa' Iscatt 



(4) 



Here p^ v ,k,aa' represents a single particle density ma- 
trix with p and v = c or v. The diagonal elements de- 
scribe the carrier distribution functions p^^.k.aa — f^ka 
of the band the wave vector k and the spin a. It 
is further noted that f ck(7 = f eka represents the elec- 
tron distribution function with a = ±~ and f vka 
I — f hka with fhka- denoting the hh distribution func- 
tion and a = ±|. The off-diagonal elements describe 
the inter spin-band polarization components (coherences) 
we defined at the end of the previous subsection with 

-iut for the inter CB-VB po- 



cv.K.aa' 



Pko 



Pk 



larization and p cc k i_i for the spin coherence. It is no- 
ticed here that for Pkaa' , the first spin index a always 
corresponds to the spin index of the electron in the CB 
(±1/2) and the second spin index a' always corresponds 
to that of the hh VB (±3/2). 

The coherent part of the equation of motion for the 
electron distribution function in the rotating wave ap- 
proximation is given by 



dt 



coh 



= dS a ilm[E»_*(t)P k(T s] + d5 a _ 4 Im[£°*(t)P fc(r _ j] 

+ 2 ^2 ^9 Im (2j P k+qaa' P k<Ja' + Pcc,k+q,-aa Pcc,k,a-a) 



gp B Blmp cc . k .a-a 



(5) 



The first two terms describe the generation rates by the 
polarized laser pulses. As mentioned before, the selec- 
tion rule requires that the optical transition can only 
happen between the conduction and the valence bands 
with the same spin direction. This selection rule is en- 
forced by the Kronecker <5-function 5 a ± i . The third term 
describes the exchange interaction correction of the ex- 
citing laser by the electron-hole attraction, thus it can 
be seen as a local field correction of the time dependent 
bare Rabi frequency dE±(t). The last term describes 
the spin flip and flop of electrons. It is noticed that if 
Impcc^cr-o- = 0, there is no spin flip and flop. Therefore 
we call Pcc,k,a-a spin coherence. It is further noticed that 
P cc fc_ii=/0* h i_i • Similarly the coherent part of the 

■ ■ 2 2 CC,n,, 2 2 

equation of motion for the hole distribution function is 
written as 

Of, 



like 



dt 



coh 



ds^ME^typ, 



■d5„*lm[E + (tyP k 



(6) 
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One notices here that differing from the electron distri- from p vv ,k,a-<7- Again this is due to the fact that 

bution function, there is no terms like gpBBlmp VVy k,a-a (Sy)±3 ± 3 = and therefore the spin of hole cannot 

in the coherent part of equation of motion for the hole be flipped. The scattering rates of f ektT and f hka for the 

distribution even if we do not neglect the contribution Coulomb scattering in the Markovian limit are given by 



Of, 



Of 



scat 



and 



2 ^ Vq2Tr5(s e k-q — £ e k + Sjk' — £jk'-q)\ (1 — f ek-qo) j eko(X ~ fjk' a') fjk' -qa' 
j—e : h,k f q(T f \ 

- fek-qa{^ ~ feka)fjk'a'0- ~ fjk'-qa') + [ R e ( Pk-qaa" Pkgg" ) 

cr" 

+ Rc(p C c,k-q,<T-aPcc,k~aa) {fjk'a' ~ fjk'-qa') + ^ Kc{Pk' a ' a" Pft - qa ' a " ) {fek-qo — /efctx) j 

2 ^ ^ Vg lirSi^Eek — q Sek ~l~ &ek' £ ek' — q)Pcc,k / .a' —a' Pcc,k' —q : — a' a' ' ^ek—qa feka) ? 
k' qa' 

— =2 V q 2 2TrS(ehk — £hk-q + £jk'-q — £ 3 k')l fhk-qa(l — fhka)(l ^ fjk'-qa')fjk'<7' 

SCat j=e,h,k'qa' (. 

- fhkaO- - fhk-qa)fjk'-qa'(l ~ fjk'a') - 2j [ Re ( P k- qa" a P kcr" a) (J 'jk' a' ~ fjk'-qa') 

a" 

+ Rc{Pk,^ q<J , a nP k > a > a >>){f hk -q a - fhka)] \ + 2 ^ Vq^SiEhk - £ h k-q + £ek'-q ~ Sek') 

J k'qa' 

Peck' — q.o' — a' Pcc,k' , — tr' a' (yfhka fhk — qcr') • 



(7) 



(8) 



One can easily prove from Eqs. (0) and (||) that ^ fe | sca t = J2 k 9 ^at° I scat = 0- This is because the Coulomb 
scattering does not change the total population of each band. 

The coherent time evolution of the interband polarization component is given by 

d i % 

jr, p kaa' = -«<W'( fc ) p w ~ ^aP-BBPk-aai + r dE_ (t) [S a l (1 - / fefe |) - p cc . k . a i}6 a ,z 

i 

+ ~ /fcfc-|) - Pccfc^-i]^'-! 

— * ^ [-Pfe+ 9 ,crcr' (1 — //ifecr' — feka) — Pk+q-aa' Pcc,k,a~a + PccM+q,a-aPk,-aa'} ■ (9) 

q 

The first term gives the free evolution of the polarization components with the detuning 

Saa'(k) = Sek + £hk - A - V q (f e k+qa + fhk+qa') (10) 

and Ao = u> — E g . Aq is the detuning of the center frequency of the light pulses with respect to the unrenormalized 
band gap. The second term in Eq. (||) describes the coupling of the optical coherence Pkaa with the forbidden optical 
coherences Pka-a and P k -aa due to the presence of the magnetic field B. This coupling makes P k a-a and P k -aa not 
small enough to be neglected in our calculation. The last term in Eq. (^|) describes again the excitonic correlations. 
The coherent time evolution of the spin coherence is given by 



d_ 



PcC,k,(T — (T 



= l -d{8 ah E-{t)P;_ a > - 5_ a iE_(t)*P k(T s) + l -d{5 a _ h E + (t)p;_ a _ 3 _ - S_ a _iE + (tyP ka _s) 

coh - 2 



+ i 



^ Vq [(/efc-j-gcr fek-\-q—a'} Pec, k, a — a Pcc.k-\-q,(T — (T {feka fek — a ) 



Pk+qacnPk-aax ~ ^k+q-aax ^kaai ] + 7i9pBB{feka ~ fek-a) 



(ii) 
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One notices from both the last terms of Eqs. (||) and (jll]) that the spin coherence causes the electrons oscillating 
between the spin-up and spin-down bands and in the mean time the unbalance between these two bands feeds back 
to the spin coherence. It is further noted from the first two terms of Eq. ( |lT| ) that pump pulse also generates the spin 
coherence. However, as this process is through the coupling to the forbidden transitions, its contribution is negligible 
compared to the effect due to magnetic field [last term of Eq. (fTl|)] as will be shown in the later sections. 
The dephasing of the interband polarization components is determined by the following scattering: 



dP, 



kacTQ 



dt 



scat 



= \ — ^ Vq2-K5(e e k~ q + £hk + £jk' — £jk'~q ~ Aq) j {Pk,aa — Pk-q,aa a ) (1 — fjk'tr')fjk'-q. 



j = c,h 

k' qa' 



^ J Pk',a'a"Pk' — q,a'a" ~t~ (/^cc,fe— q,a— aPk,— tjtJo fek — qaPkaao fhka^Pk — qaa^i^fjk' a' fjk' — qa' ) J 1 



Vq2nS(e e k- q + £hk + £ e k' — Sek'-q — &o)(Pkaao ~ Pk-qoo ) Pcc,k' ,a' -a' Peck 1 

k' qa' 

Pkaao 



{k k — q, k k' — q} 



T-2 



(12) 



Here T2 is introduced phenomenologically to describe additional slower scattering process such as carrier-phonon 
scattering, {k <-> k — q,k' «-> k' — q} in Eq. ( |l2] ) stands for the same terms as in the previous {} but with the 
interchanges k k — q and k' «-> k' — q. We point out here that all the scattering terms in Eqs. (0), (ra) and ([jj]) are 
the contributions from the Coulomb scattering. We did not include the EHSE scatterings as they are much weaker 
in comparison with the Coulomb scatterings. 

The Coulomb scattering contribution to the scattering term of the spin coherence is: 



dpcc,ka — a 



dt 



Coul 



scat 



< - ^ V q 2 2irS(e e k-q - S e k + £jk> ~ £ 3 k'-q)\{fek-qa 

I j = i,h 

qk' a' 



— qaPcc,k,a—o 



Pcc,k — q,a—afek—a 



^ t Pk-qoo" Pk— a-a" ) (fjk' a' fjk' -qa') + Pcc,k,a-a (1 fjk' a') fjk' -qa' ^ ' Pk' a' a" Pk' -qa' a 



ePcc. k' ,a' —a' Pcc,k' —q,—a' 



Pcc,k — q,a — a 



fjk'a'i^ fjk'-qa') ^ ] Pk' a' a" Pk' -qa' a" 



Sj=ePcc,k',a'—a'Pcc,k'-q,-a'a' — {fc<-^fc — q, k ^ k q} . 



(13) 



One can prove from Eq. ([13]) analytically that 
a Coul 

ClPcc.ko-a 
k gca t 







(14) 



This can be easily seen as the second half of Eq. fll3| ) 
comes from the first half by interchanging k <-> k — q and 
k' <-> k! — q. In the sum over k of Eq. (|l4|), one may per- 
form the following variable transformations: k — > — k + q 
and k' — > — fc' + g to the second half of Eq. (|l3l ) which 
make the second half just the same as the first one but 
with opposite sign. Eq. ([l4[) indicates that the Coulomb 
scattering does not contribute to the spin dephasing. 

In order to study the dephasing of the spin coherence, 
we pick up the contributions from EHSE scattering. The 
"direct" EHSE contribution is given by 

dEHSE 



dp, 



cc.kcr—a 



dt 



= { ~ Jq ^ Uq27r5(s e k-q ~ S ek ~ S h k' + Zhk'-q) 
k' qa' 

^ [(^ ~~ feka fek — a)Pcc,k — q,a — afhk' — qa'{X fhk'a') 
Pec k, a— a {fek — qa fek — q — a ) fhk' — qa' ( 1 ~ fhk' a' ) 

+ {k^ k-q, k' <-> fc'-g} . (15) 
The "exchange" EHSE contribution can be written as 
eEHSE 



dpec, 



ka — a 



dt 



scat 



{"jg ^ ^q 2 27T^(e efc / + £hk'-q — Sek ~ £}lk-q) 



k' qa' 



scat 



^ [(^ fek'a fek'—a)Pcc,k,a — afhk — qa'(^ fhk' — qa') 
~t~ Pcc,k' ,a — a (feka ~l~ fek — a ) fhk — qa' (1 fhk'—qa' ) } 

+ {fc^/c'}. (16) 
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It is noted that in Eqs. ( |l5| ) and (|T^), the second 
half of each equation shares the same sign as the 

dEHSE 



first half. 



Therefore J2k 



at 



scat 



^ and 



at 



cEHSE 

^ 0. This tells us that EHSE con- 
scat 

tributes to the s pin dcphasing. One further finds by 
comparing Eqs. (|l^) and (|l6|) that besides the above 
mentioned U' q <^U q , the energy phase space of the "ex- 
change" EHSE, which is imposed by the energy conserva- 
tion, is quite limited compared to that of "direct" EHSE. 
All these indicate that the contribution of the "exchange" 
EHSE is negligible in comparison with that of the "di- 
rect" EHSE. We have further proved that there is no 
contribution to the scattering term from the combina- 
tion of Coulomb scattering and EHSE as U q V q or U' q V q . 
Therefore 



dp c 



,kcr—c 



dt 



dpc 



,ka—c 



scat 



dt 

dp C c,ko 



Cou 



dp c 



.k(7—c 



scat 



dt 



dEHSE 



scat 



dt 



eEHSE 

scat 



(17) 



Equations (|)-@ comprise the complete set of kinetic 
equations. It is noted that we only include EHSE in the 
scattering term of the spin coherence. Its contribution to 
the optical coherences and electron (hole) distributions 
are neglected as it's much smaller than the contribution 
from the Coulomb scattering.^ Moreover, the electron- 
hole recombination is not included in our model as the 
time scale for such effect is at least one order of magni- 
tude longer than the time scale of dephasing. 



C. Faraday rotation angle 

The FR angle can be calculated for two degenerate 
Gaussian pulses with variable delay time r. The first 
pulse (pump) is circular polarized, eg. -Epump = E_(t), 
and travels in the direction ki . The second pulse (probe) 
is linear polarized and is much weaker than the first one, 



X [E°_(t 



E° + (t 



t)] with ^< 1, The probe pulse 



travels in the k2 direction, .— 
The FR angle is defined asci? 



^f(^probJ^) 

-Pk^m^Jt-r^dt, (18) 

with Pkaa standing for the optical transition in the prob 
direction, i.e. k2 direction. C is a constant. 

For the delay time r is shorter than the optical de- 
phasing time, one has to project the optical transition 
Pkaa to k2 direction. One may use an adiabatic pro- 
jection technique described in detail in Ref. E9. This 



technique is suitable for optically thin crystals, whepa 
the spatial dependence can be treated adiabatically.E£l 
To do so, one replaces the single-pulse envelope func- 
tion in Eqs. (||) and (g) by two delayed pulses E°_(t) = 
E°_(ty* + E° woh _(t - r) and E° + (t) = E^ Jt - r) 

with the relative phase tp = (ki — k2) • x resulting from 
the different propagation directions. The projection tech- 
nique is used with respect to this phase. However, when 
delay time r is much longer than the optical dephasing, 
the optical transition Pkaa induced by the pump pulse 
has already decayed to zero and therefore one may per- 
form the calculation with tp = 0. 

It is interesting to see from Eq. ( |l8| ) that although the 
spin coherence is determined by p cck i_i, it does not 
appear directly in the final equation of the FR angle. 
Instead, it affects the FR angle through the optical tran- 
sitions Pkaa ■ For delay time t much longer than the opti- 
cal dephasing time, the optical coherences Pkaa induced 
directly by the pump pulse together with the forbidden 
optical coherences Pka-a have already been destroyed to 
zero. However, the spin coherence induced by the same 
pump pulse remains and makes the electrons oscillate 
between the spin- up and down bands. This unbalance 
of population strongly affects the optical transitions in- 
duced later by the probe pulse around t and gives rise to 
the time evolution of Faraday rotation angle. 



III. NUMERICAL RESULTS 



We perform a numerical study of the Bloch equations 
to study the spin coherence of optically excited electrons 
in an undoped insulating ZnSe/Zni_. E Cd ;r Se quantum 
well. As the main interest of the present paper is fo- 
cused on the mechanism of the spin dephasing, we will 
not perform a pump-probe computation to calculate the 
FR angle Eq. (|l^) as it requires extensive CPU time-,and 
also as it has been calculated by Linder and Shamlla for 
the same set of Bloch equations but with relaxation time 
approximation for all the dephasing. Instead, we will 
only apply a single pump pulse and calculate the time 
evolutions of both the optical and spin coherences to- 
gether with the electron and hh distributions after that 
pulse under the carrier-carrier scattering. The dephasing 
of the spin coherence is well defined by the incoherently- 
summed spin coherence, p{t) = J2k \Pcc,k,±-±{t)\i as we ^ 
as the optical dephasing is described by the incoherently 
summed polarization, P a a{t) — J2k \Pk,<r<Mi)\- The later 
was first introduced by Kuhn and Rossi.E! It is under- 
stood that both true dissipation and the interference of 
many k states may contribute to the decay. The incoher- 
ent summation is therefore used to isolate the irreversible 
decay from the decay caused by interference. From p(t) 
and Paa(t) one gets the true irreversible dephasing of 
spin and optical coherences respectively. 



G 




FIG. 1. Electron distributions vs. time t and electron en- 
ergy E e for the spin 1/2 (Fig. 1(a)) and -1/2 (Fig. 1(b)) and 
hh distribution vs. t and hh energy Eh (Fig. 1(c)). 

The material parameters in our calculation are taken 
from the experimental data with effective mass m e = 



0.152mo for electronl23 and rrih = 6m e for hhJ2a Exci- 
ton Rydberg En = 19 meV and the g-factor is take 
as 1.3. This abactor is larger than what reported in 
the experimentEIJ (1.1) by measuring the beat frequency 
of the FR angle. The reason will be shown clearly in 
later subsections. is taken as 10 ps for all the cal- 
culations. We choose a left circular polarized Gaussian 
pulse E°_(t) with the width St = 100 fs and detuning 
A = 4.275 meV. The total density excited by this pulse 
is 3 x 10 11 cm _ ™which is the same as what reported in 
the experiment O E^_(t) = 0. Under such a pulse, one 
can see immediately from the kinetic equations that only 
hh with spin | can be excited and fhk-^ = 0- Moreover, 
P h i 3 = P h i 3 = 0. Therefore, we only need to deter- 

K 2~2 ^2 2 

mine the electron distribution functions f e ka (cr = ±1/2), 
the hh distribution functions f^i, interband polariza- 
tions (optical coherences) Pfci| and P fe _i 3, and spin co- 
herences p cc i_i. We will solve the Bloch equations with 
only Coulomb scattering and with both Coulomb scat- 
tering and EHSE scatterings respectively. 
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FIG. 2. Absolute amounts of optical coherence |P fc i3(t)| 
[Fig. 2(a)] and spin coherence \p cc k i _ i (t)| [Fig. 2(b)] as 
functions of t and electron energy E e . Note the time scale 
of optical coherence is much shorter than the spin coherence. 



A. Coulomb scattering 

We first discuss the Boltzmann kinetics under the 
Coulomb scattering with the statically screened instan- 
taneous potential approximation: 



V a 



q e a (q + K) 
The inverse screening length k is expressed asM ! 1 



(19) 



= —Y[U=^{t) + -±f hk ^{t)\ . (20) 
e m e 



It is noted that this static screening model, although sim- 
plifying numerical calculation significantly, overlooks the 
plasmon-mode contributions and overestimates the elec- 
tron screening. 
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FIG. 3. Total densities of each spin band N e kcr{t) for elec- 
tron and N hk 3,(t) for hh (solid curves) together with the in- 
coherently summed polarization Pi 3 and spin coherence p(t) 
(dashed curves) are plotted against time t for B — 4 T. Note 
the scale of the coherences is on the right side of the figure. 

We first show in Fig. 1 the electron and hole distribu- 
tion functions f e ka{t) (f = ±1/2) and fh k i(t) versus t 
and electron energy s ek (for electron distributions) or hh 
energy e^k (for hole distribution) for B = 4 T. In the 
initial time one observes a small peak around 10 meV for 
f ek i in Fig. 1(a). This peak is the effect of the pump 

pulse and the strong Hartree correction Eq. (10). The 
similar peak can also be observed for the hh distribution 



function in Fig. 1(c). At later times the carriers relax to 
the low energy states and f ek i reaches its first highest 
peak at 1.4 ps. Around 7 ps, the distribution function 
of spin-up band reaches the valley of zero values and the 
spin-down band, in the mean time, arrives at its peak. 
This indicates that electrons in the spin-up band evolve 
into the spin-down band. After that electrons start to 
move back to the spin- up band and at about 14 ps f ek i 
reaches its second highest peak and / e fc-i reaches its 
valley. The reason that the second highest peak is a lit- 
tle higher than the first one is because there are more 
electrons at higher energy states at initial times due to 
the position of the pump pulse. This oscillation keeps 
on without any decay. The distribution for hh relaxes 
into the Fermi-like distribution after a few picoseconds 
and remains unchanged. In Fig. 2 the absolute value of 
the optical coherence, \P k ia(t)|, and the absolute value 
of the spin coherence, \p cc k i_i (t)\, are plotted as func- 
tions of t and electron energy. It is seen from Fig. 2(a) 
that the optical coherence decays very quickly and within 
first few picoseconds it has already totally disappeared. 
Nevertheless, the spin coherence does not decay at all. 
The second peak is of the same height as the first one. 

The incoherently summed polarization, Pis(t), and 
the incoherently summed spin coherence, p(t), are plot- 
ted as dashed curves in Fig. 3. The total densities of 
each spin band N ea (t) — ^2 k f e ka(t) for electron and 
N h z(t) = ^feAfef(^) f° r un are a l so plotted as solid 
curves in the same figure. It is seen from the figure that 
the optical coherence injected by the pump pulse is about 
3 times larger than the spin coherence. However, this 
coherence is strongly dephased by the Coulomb scatter- 
ing and is totally gone within the first few picoseconds. 
It is further shown from the figure that the spin coher- 
ence exhibits beating which does not decay at all. The 
electron densities of spin-up and down bands oscillate 
between zero and the total excitation. These results con- 
firm that pure Coulomb scattering does not contribute to 
the spin dephasing. We further find that the frequency of 
the oscillation is mainly determined by the Zeeman split 
gj-isB, but red-shifted by the Hartree- Fock terms in Eq. 
(|ll|). The reduced ^-factor resulting from the oscillation 
frequency is 1.25. It is interesting to see from the figure 
that the maximum value of |iV e i —N e _i \ occurs when the 
(incoherently summed) spin coherence is zero. The for- 
bidden optical coherence is about 30 times smaller than 
the optical coherence and decays similarly as the optical 
coherence. It plays an insignificant role in this problem. 



B. EHSE 

Now we include the contribution of EHSE [Eqs. Jl5| ) 
and ([l6])] . As said before that the contribution from "ex- 
change" EHSE is negligible, we therefore only include the 
"direct" one [Eq. Jig)] here. We are lacking information 
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of the matrix elements U q which requires a detailed band 
structure calculation. For the sake of simplicity and also 
for comparison with the effect of Coulomb scattering, in 
this study we assume it is the same as V q but with a 
phenomenological prefactor 4\^F/3. We take T = 0.015 
so that in the scattering term, Eq. ([l5]), the matrix ele- 
ment of EHSE is two orders of magnitude smaller than 
that of the Coulomb scattering. This number is taken 
so that the sain dephasing time is in agreement with the 
experiment O We have also performed numerical calcu- 
lation by taking only the "exchange" EHSE and found in 
order to get the same spin dephasing, the prefactor has 
to be 0.1 which is one order of magnitude larger than 
the present case. This big prefactor is understood due 
to the effect of the limitation of the energy phase space 
discussed before. 




FIG. 4. Electron distributions vs. time t and electron en- 
ergy E e for the spin 1/2 (Fig. 4(a)) and -1/2 (Fig. 4(b)) for 
B = 4 T. Effects of EHSE are included. 

The electron distribution functions f e ka{t) {p = ±1/2) 
versus t and electron energy e e k are plotted in Fig. 4 for 
B = 4 T. The hh distribution function /^a(t) versus 



t and hh energy etk remains unchanged from Fig. 1(c) 
after the inclusion of EHSE. In the initial time one ob- 
serves again a small peak around 10 meV for f ek i in 
Fig. 4(a) due to the effect of pump detuning. A similar 
peak can also be observed for the hh distribution func- 
tion as in Fig. 1(c). Again one observes that the carriers 
relax at later times to the low energy states and all the 
distributions show the Fermi-like distributions. Differ- 
ing from the case with only Coulomb scattering, there 
are only small oscillations of electrons between spin 1/2 
and —1/2 bands in the later time. The hh distribution 
remains unchanged after the first few picoseconds as be- 
fore. The absolute value of optical coherence, |P fc ia(t)|, 
versus t and electron energy is also unchanged from Fig. 
2(a) after inclusion of EHSE. Therefore, as before, the 
optical coherence decays very quickly and within a few 
picoseconds it has totally disappeared. However, EHSE 
makes big change to the absolute value of the spin co- 
herence as can be seen in Fig. 5 where \p cc k is 
plotted as a function of t and electron energy. It is seen 
that the spin coherence lasts much longer and one can 
see a smaller second peak around 10.5 ps and a much 
smaller third peak around 18.4 ps. Comparing with Fig. 
2(b), one can see the effect of strong spin-dephasing by 
EHSE as the second peak is much lower than the first 
one. For t > 20 ps, the spin coherence has almost gone. 
One more point needs to be addressed here is, one can 
see from Fig. 5 that there is a very small bump around 
10 meV at initial time. A similar bump can also be seen 
in Fig. 2(b). These bumps are the effects of the pump 
pulse described before after Eq. (pd|). One can see that 
they are much smaller compared to the effects caused by 
the magnetic field. This also justifies the assumption we 
made before that hh-hh spin coherence can be neglected 
in this investigation. 




FIG. 5. Absolute amounts of spin coherence \p cc k 11 (i)| 
as functions of t and electron energy E e . Effects of EHSE are 
included. B=4T. 







The incoherently summed polarization, Pis(t), and 
the incoherently summed spin coherence, p{t), are plot- 
ted as dashed curves in Fig. 6. In order to have the effect 
of spin coherence more pronounced, we plot in the same 
figure also 4.5 x p(t) as dash-dotted curve. The total 
densities of each spin band N ea {t) and N h z{t) are also 
plotted as solid curves in Fig. 6. It is seen from the fig- 
ure that the optical coherence injected by the pump pulse 
is about 4.5 times larger than the spin coherence. How- 
ever, this coherence is strongly destroyed by the Coulomb 
scattering within the first few picoseconds. This is part 
of the reason of the fast decay in the ER angle for the 
first few picoseconds in the experiment O As seen in Eq. 
(|l8|), the FR angle is proportional to the "total" opti- 
cal transitions. For the first few picoseconds, the optical 
transitions are induced by both the pump pulse and the 
much weaker probe pulse, and the fast decay of the ini- 
tial pump pulse-induced optical transition gives rise to 
the strong decay of the FR angle. 
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FIG. 6. Total densities of each spin band N e k a {t) for elec- 
tron and N hk 3.(t) for hh (solid curves) together with the in- 
coherently summed polarization Pi a and spin coherence p(t) 
(dashed curves) are plotted against time t for B — 4 T. The 
dash-dotted curve is 4.5p(i). Effects of EHSE are included. 
Note the scale of coherences is on the right side of the figure. 



It is further shown from the figure that the electron 
densities of spin up and down bands oscillate with the 
amplitude of the oscillation decaying. The circularly po- 
larized pump pulse first pumps electrons into the spin-up 
CB from the spin-up hh band. Therefore electrons first 
occupy the spin-up CB band and leave behind hhs in the 
spin- up VB. Therefore N e i and N h 3 fast rise to 3 x 10 11 
cm~ 2 within the time scale of the pump pulse. Then 
due to the magnetic field, electrons in the spin-up band 
start to go to the spin-down band. This makes N e i de- 
creases and N~_ i rises. In the mean time, the unbalance 



in populations also serves as pump field to the spin coher- 
ence. After t is around 5 ps, the electron population in 
the spin-down band surpasses that in the spin-up band. 
Without dephasing, this oscillation may keep going on as 
shown in the previous section. However, the spin dephas- 
ing makes these two populations finally merge. From the 
figure, one can see that for t > 20 ps, the difference is al- 
ready negligible compared to the oscillations before and 
p also decays to zero. This oscillation is also shown in 
the experiment through the FR angle as beatings for the 
same magnetic fieldO All these results confirm what we 
proposed in Sec. II, that the spin dephasing is caused by 
EHSE and the electron spin coherence is represented by 
p. Moreover, one also finds that besides the effect of spin 
dephasing, EHSE also changes the oscillation frequency 
and hence the beating frequency of the FR angle. By 
comparing the period of the oscillations in Figs. 3 and 
6, one finds the period increases by about 1.34 ps in the 
later case. This means that the EHSE further red-shifts 
the spin sphtting. Therefore, the g factor reported in the 
experiment^ by measuring the frequency of the beating 
of the FR angle is the effective g factor. It is seen from 
the figure that the period of the oscillation is T = 15.58 
ps. As 2ir/T = g e QpBB, the effective g-factor g e g- is 
therefore 1.14, in agreement with the what reported in 
the experiment as 1.1. EU Again from the Figure one ob- 
serves that the maximum value of |iV e i — iV e _i | occurs 
when the (incoherently summed) spin coherence is zero. 
Our calculation shows again that the forbidden optical 
coherence is insignificant in this problem. 
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FIG. 7. Total densities of each spin band N e k a (t) for elec- 
tron and N hk a(t) for hh (solid curves) together with the in- 
coherently summed polarization Pi 3 and spin coherence p(t) 
(dashed curves) are plotted against time t for B — 2 T. Ef- 
fects of EHSE are included. Note the scale of coherences is 
on the right side of the figure. 



In order to compare with B — 2 T case, we plot in 
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Fig. 7 the incoherently summed polarization Pi 3 (t) and 
spin coherence p(t) as function of time t, together with 
the total densities of each spin band N ea (t) and N h a(t). 
One finds the spin coherence decays at the same rate as 
in the B = 4 T case but without any beating. So does 
the electrons in the spin bands. This confirms the finding 
in the experiment that there is no beating in the Faraday 
rotation angle.EJ 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we have performed theoretical studies of 
the kinetics of the spin coherence of optically excited elec- 
trons in an undoped insulating ZnSe/Zni-^Cd^Se quan- 
tum well under moderate magnetic fields in the Voigt 
configuration. Based on a two spin-band model in both 
the CB and the VB, we build the kinetic equations com- 
bined with intra- and interband Coulomb scattering and 
interband EHSE in the Boltzmann limit. We include all 
the coherences induced directly by the laser pulse — opti- 
cal transitions — and indirectly through the effect of the 
magnetic field — electron-electron spin coherence and for- 
bidden optical coherence — for the electron and hh in our 
model. The hh-hh spin coherence p vv k 3_3 is neglected 
in our present investigation because it can not be induced 
by the effect of the magnetic field, but only through the 
pump pulse coupled to the forbidden transitions and is 
therefore much smaller than the other coherences. We 
separate the spin coherence from the well known opti- 
cal coherences and study the effects of Coulomb scat- 
tering and EHSE on all the coherences. We find that 
the Coulomb scattering makes strong dephasing of the 
optical coherence and forbidden optical coherence. How- 
ever, it does not contribute to the spin dephasing at all. 
EHSE is the main mechanism leads to the spin deco- 
herence. We numerically solve the kinetic equations for 
two different magnetic fields. We find that the beating 
in the Faraday rotation angle is basically determined by 
the electron Zeeman splitting gpsB, however, with a red 
shift from the Coulomb Hartree-Fock contribution and 
EHSE effect. The forbidden optical coherence is found 
of marginal importance in this problem. The matrix el- 
ement of the Coulomb scattering is taken as statically 
screened Coulomb potential and the matrix element of 
EHSE is assumed the same as the Coulomb scattering 
with a phenomenological prefactor. By fitting this-pref- 
actor with the dephasing time in the experiment JIj our 
theory can well explain the experiment data for two dif- 
ferent magnetic fields. A first principal investigation 
of the EHSE scattering matrix clement is definitely im- 
portant fo r a m ore thorough understanding of the spin 
dephasing.EIE3 

For the n-doped material, things are quite different 
from the undoped case discussed in this paper. Due to 
the presence of large numbers of doped electrons (about 
one order of magnitude larger than the optically excited 



carriers) in the CB, the lifetime of holes is therefore quite 
short and is measured smaller than 50 ps compared to the 
lower limit of 100 ps for the insulating sample.EJ The ex- 
periments found that the spin dephasing for the doped 
sample is three orders of magnitude longer than the un- 
doped sample, but with a fast dephasing (losing about 
50% coherence) within the first 50 ps. The fast dephas- 
ing in the initial times can be well understood by the 
EHSE discussed above, but modified with the fast de- 
creasing of hh population. Nevertheless, the mechanism 
of the spin dephasing discussed above cannot be applied 
to the doped case after 50 ps as shown from Eq. ( |l5| ) 
that the scattering of EHSE is proportional to the hole 
distribution and therefore its contribution to the spin de- 
phasing decreases with the recombination of electron and 
hole. However, for the n-doped case, it is also possible 
for EHSE to contribute to the spin dephasing. As the 
timescale of the spin dephasing for the n-doped sample is 
much longer than the insulating sample, the coupling be- 
tween the hh and the light hole can not be neglected and 
the spin of hh also experiences the precession. Eil There- 
fore, the hh-hh spin coherence p vvk ^_^, which is in- 
significant for the undoped case, may play an important 
role in the doped sample. Its contribution to the scatter- 
ing term to the electron spin coherence is in the similar 
form as Eq. (|l5|) but all the hole distribution parts are 
replaced by terms composed of hh-hh spin coherences. 
In the absence of the hole distribution, it becomes the 
leading mechanism from the contribution of EHSE to 
the spin dephasing. Physically this contribution is the 
spin exchange between electrons and virtual holes. This 
mechanism, together with other spin dephasing mecKa-, 
nisms due to the band mixing for the n-doped case,EHI 
is still under investigation. A corresponding extension 
of the present theory for the n-doped material will be 
published in a separate paper. 
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This many body interaction Hamiltonian can be de- 
rived from the interaction potential - 



[~2 



8 -r5(r)] 



+ 



JzS z . The first term 



r=|ri-r 2 | 

describes the Coulomb potential between two electrons. 
One can use a plane-wave approximation to get-jthe first 
term in Eq. (||). The second term is the EHSEEI with J 
representing the spin of the electron in the hh VB and 
S denoting the spin of the electron in the CB. For the 2D 
case, a plane-wave approximation cannot be applied to this 
term and one has to use the wave function of the material 



to calculate the matrix elements [/<- 



and Uq for the second 



and third terms in Eq. (|3|). 

It is noted that it is also true that the contribution from 
EHSE to the spin-coherence scattering term is also much 
smaller than that from Coulomb scattering. However the 
EHSE 



fact that y"] fr 



scat 



7^ makes it contribute to 



the long time decay of the spin coherence whereas the fact 
EHSE 



at 



scat 



= makes its contribution to the 



that y\ 

scattering of electron (hole) distribution negligible. Its con- 
tribution to the scattering of optical coherence is obviously 
negligible as the later vanishes within first few picoseconds. 
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